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We study a closed model of a universe filled with viscous fluid and quintessence matter com- 
ponents. The dynamical equations imply that the universe might look like an accelerated flat 
Friedmann-Robertson- Walker (FRW) universe at low redshift. We consider here dissipative pro- 
cesses which obey a causal thermodynamics. Here, we account for the entropy production via 
causal dissipative inflation. 

PACS numbers: 04.20Jb,98.80.Cq 



INTRODUCTION 



Recent observational evidences suggest that the measured matter density of baryonic and dark matter is significantly 
less than one, i.e. its critical value. This implies that either the universe is open or that there are some other 
matter components which makes fltotai ~ 1- Combined measurements of CMB temperature fluctuations and of the 
distribution of galaxies on large scales strongly suggest the possibility of a flat universe Q, , which is consistent 
with the standard inflationary prediction 3] . Recent measurements of a la type distant supernova (SNela) [ill, at 
^ ■ redshift z ~ 1, indicate that the expansion of the present universe is accelerated. 

An initial interpretation was to consider that in the universe there exists an important matter component that, in 
its most simple description, has the characteristic of the cosmological constant A, i.e. a vacuum energy density which 
contributes to a large component of negative pressure, and thus accelerates rather than decelerates the expansion of 
the universe. An alternative interpretation is to consider quintessence (or dark energy), in the form of a scalar field 
with a self-interacting potential. In the literature various examples of such Q a component have been considered 
starting from the pioneering paper which also considered the case of cosmic strings. The non-minimally coupled 
scalar fields were considered in Ref. Q and the scalar-tensor theories in Ref. @|. The current status about dark 
energy can be found in Ref. @. 

On the other hand, in recent years important attention has been received by cosmological models which consider 
bulk viscosity. For instance, it was shown that the introduction of this kind of viscosity into cosmological models can 
avoid the big bang sin gula rities [ToL ITU , and any contribution from particle production may be modeled as an 



i gula 

effective bulk viscosity [f| ■ The bulk viscosity typically arises in mixtures of different species (as in a radiative fluid) 
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or of the same (but with different energies) fluids. The dissipation due to bulk viscosity converts kinetic energy of the 
particles into heat, and thus one expects it to reduce the effective pressure in an expanding fluid. This fact may play 
a crucial role in the inflationary era of the universe since it is interesting to know whether dissipative effects could be 
strong enough to make a large negative effective pressure leading to inflation. 

Many, probably most, of the inflationary cosmological models considered are of a flat FRW type, since they are 
spatially isotropic and homogeneous. This implies that velocity gradients causing shear viscosity, and temperature 
gradients leading to heat transport, are absent. Thus any dissipation in an exact FRW universe is scalar, and therefore 
may be modeled as a bulk viscosity within a thermodynamical approach. These scalar dissipative processes may be 
treated in cosmology via the theory of general relativity, considering the bulk viscosity [lil ITlj , which is compatible 
with the homogeneity and isotropy assumptions for the universe. These dissipative processes may play an important 
role in the early universe, especially before nucleosynthesis ^(|. Then we shall use the causal thermodynamical 
theory for processes not in equilibrium |17| . The stable and causal thermodynamics of Israel and Stewart replaces 
satisfactorily the unstable and non-causal theory of Eckart and Landau and Lifshitz [l^ . 

It is well know that the inflationary universe model is necessary for solving most of the cosmological puzzles, in the 
way that it explains how a large class of initial states can evolve into a unique final state that is consistent with our 
observed universe. In most of inflationary FRW cosmological models not only the scale factor a(t) presents a positive 
acceleration a > 0, but also it is assumed from the very beginning that the geometry of the universe is completely 



In the light of this approach an interesting question to ask is whether this flatness may be due to a sort of 
compensation among different components that enter into the dynamical equations. In this respect, one of our goals 
is to address this question in a frame of a simple model. In the literature we find some descriptions along these 
lines. For instance, a closed model has been studied with an important matter component whose equation of state 
is given by p = — p/3. There, the universe expands at a constant speed j2jj. Other authors, while using the same 
astronomically observed properties for the universe, have added a nonrelativistic matter density in which the total 
matter density Qo is less than one, thus describing an open universe [2l]]. Also, a flat decelerating universe model 
has been simulated 22]. The common fact to all of these models is that, even though the starting geometry were 
other than that corresponding to the critical geometry, i.e. a flat geometry; all these scenarios are, at low redshift, 
indistinguishable from a flat geometry. 

In this paper we wish to consider a closed universe model composed of two matter components: one related to 
a viscous matter and the other to quintessence matter. The geometry, together with these matter components, 
confabulates in such a special way that it gives rise to a flat accelerating universe scenario. In this regard, we use 
the standard anzats (see Eqs. I|16|) . (|I4() and l|24|) below) which introduce five parameters in our model, where there 
exists a range of these parameters that gives an appropriate value for the entropy of the universe. The value obtained 
for the entropy in our model agrees with the expected value obtained from usual inflationary re-hea ting . In this way, 
as we will see, our model accounts for the generally accepted entropy production, via warm inflation[23j. In this kind 
of model the radiation is continuously produced by the decay of the inflaton scalar field. In this way this field should 
be coupled to ordinary matter and, contrary to usual inflation, primordial density fluctuations are originated much 
more from thermal fluctuations than from quantum fluctuations of the inflaton 24]. At the end, the parameters of 
the resulting model could be estimated by using astronomical data. 

The outline of the present paper is as follows: In Sec. II we briefly review the system of gravitational field equations 
and the transport equation for bulk viscosity. In Sec. Ill and IV a power law FRW and de Sitter FRW causal 
cosmologies are considered. In Sec. V some conclusions are given. 

II. THE SYSTEM OF GRAVITATIONAL FIELD EQUATIONS AND THE TRANSPORT EQUATION 

FOR BULK VISCOSITY 

Now we shall obtain the field equations for the considered cosmological viscous model. The energy-momentum 
tensor of the fluid with a bulk viscosity is given by 



where qm and pm are the thermodynamical density and the pressure of the fluid, while it is the bulk viscous pressure, 
and i/' = Sq. 

The contribution to the energy-momentum tensor related to the classical scalar field 0, which is minimally coupled 
to gravity, becomes given by: 



flat. 



Tfiu = (pM +PM + rfVpV,, + (p M + n)g, 



(1) 




(2) 
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where V((f>) is the scalar potential. 
Now we consider the FRW metric 



ds 2 = -dt 2 + a 2 (t)[ ^ 5 + r 2 (d9 2 + sin 2 9 d<j) 2 ) I . ( :! ) 

1 — fcr 2 



From Eq. and © we have the gravitational field equations: 



3H 2 + 3^ =k (\tf + V + p M ) (4) 



2ij + 3ff 2 + 4r = K (-^ 2 + ^-p M -^j (5) 

where H = and k is the curvature parameter and n — 8nG, with G the Newtonian gravitational constant. 
The fluid conservation equation is given by 

p u + 3H(p M + p M + 7r) = (6) 

and the equation for the scalar field by 

dV 

+ 3if^+— =0. (7) 
The energy density associated with the scalar field is given by 

P* = \& + V{4>) (8) 



and the pressure by 

p, = \4> 2 - V(4). (9) 

The scalar field has a state equation defined by p . = wp^ and satisfies the conservation equation 

p^+3H(p^+ P J=0. (10) 

^From Eq. @ and the weak energy condition p^ > we conclude that, for desired flatness simulation, i.e. Eq. Q 
becoming the well known Friedmann Equation for a flat FRW universe 

3H 2 = k Pm , (11) 

we must take k = 1 in the field equations. Thus 

*P* = 3/« 2 - (12) 

^From here we have that Kp^ = —3a /a 3 = —6H/a 2 and, substituting into Eq. (|10fl . obtain that w = — |. Notice that 
this result is independent of the explicit expression for the scale factor as a function of t. 

At this point we should notice that our model is intrinsically closed, therefore we should stress that our model is 
completely different with the flat standard Einstein model which characterizes with VLm = 1 and is ruled out by the 
currents astronomical observations. In fact, the curvature term (which is intrinsically geometric) cannot be eliminated 
by a simple redefinition of the scalar field tf>, since the curvature is a geometrical property, which follows directly from 
the metric tensor and which enters into the FRW line elements. These two models become indistinguishable only fol- 
low redshifts and are similar to the cases studied in Refs. [2ll l22l |2{| . 

Another point that we would like to stress here is that almost all physical quantities are obtained in terms of the 
scale factor and its derivatives. From Eq.^]we get p^ as a function of a. Since w = —1/3, then <\? = V((f>) — 2/na 2 . 

Also, from Eq.|S]we get that ktt = —23H 2 — 3jH 2 , where we have used the state of equation pu = (7 — l)p M for the 
matter component and Eq. 1111 As we will see and in order to obtain an explicit expression for the scalar potential V 
as a function of the scalar field (j>, we need the temporal dependence of the scale factor. In this respect we will take a 
power law and an exponential inflationary solutions. 
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^From now on we shall use the simplest case of scalar dissipation due to bulk viscosity £, modeling cosmological 
dissipation from the consideration of noncquilibrium thermodynamics via Israel-Stewart theory. The bulk viscous 
pressure ir is given by the transport equation (linear in n): 

T7r + 7r = -3Jr£-iTTr^3Zr+£-|-!^, (13) 

where r is the relaxation time (which removes the problem of infinite propagation speeds), £ the coefficient of bulk 
viscosity and T is the temperature of the fluid. In the non-causal formulation r = and then Eq. (|13|) has a simple 
form 7r = —3H£. 

Following 0, |2(| we shall take the different thermodynamic quantities to be simple power functions of the density 

£ = ap M m , T = fip M r and r = -L = ap,/"" 1 (14) 

Pm 

with a,n,m and r positive constants. The expression for r is used as a simple procedure to ensure that the speed of 
viscous pulses does not exceed the speed of light. For an expanding cosmological model the constant m should be 
positive and 

t>H- x (15) 
in order to have a physical behaviour for the coefficient of viscosity £ and the relaxation time r [Til [27J ■ 

III. POWER LAW FRW CAUSAL INFLATIONARY COSMOLOGY 

Now we are interested in power law inflationary models; thus we shall take the scale factor in the form 

a(t) = a (£\ , (16) 
where clq and n are constant parameters, with n > 1. Then, using Eq. !)!()[) we obtain 



3 Ao x2 



ao V t 

and 

ao 2 \t 



kP^ — It) ( 18 ) 



^ = -^(7) (19) 



for the Hubble constant, mass density and pressure of the field <j> respectively. 
In this case the solution for the scalar field and its potential are given by 

m - ^0 ( j) (20) 

and 

V{4>) = V (J^j " _I , (21) 

respectively, where <p = <j>{t ) = and Vb = 

We shall consider for the fluid a barotropic equation of state 

P„ = (7-1)P„, (22) 
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where < 7 < 2. 

Then for the power law expansion and from eqs. ©, (|13fl and l|14f> (with the standard thermodynamic relation 
for the temperature of a barotropic fluid r = (7 — l)/7) we obtain 

Pu + (3n + 1)^ + -p l - m PM + ^ P l7 m 
t a at 

_ 9n2 f 1 _l)^L_(^_l)^l = 0. (23) 
V 2/ t 2 2 7 p M v ; 

Note that this equation may also be written in terms of H using Eq. pill . 
Case m =/= 1: For solving Eq. (|23|1 we shall take the following anzats|2?j| : 

'»-*(£) ■ (24) 

For physical sense we must have m < 1. Inserting (|24|l into Eq. (|23[) we obtain, for m ^ 1, that the constant is 
given by 



/ (27)- 1 - 3n(l - m) - 9n 2 (l - 7 /2)(l - to) 2 \ 1/(1_m) 
^ (l-3n 7 (l-m)) J 

Later on we shall assume that the parameters 7, n, m and are given; then we can find the parameter a. The 
pressure p is obtained by using the equation of state 12211 . 

Now from Eq. H22fl and Eq. (B) we obtain for the bulk viscous pressure 



*=-%-1P» ( 26 ) 



and substituting (|24|l and p M into (|26[) we obtain 



1 \ n /t- 1/(1 - 



3n(l - m)y \t y 

^From here we see that 7r < if 7 > l/(3n(l — to)), vr > if 7 < l/(3n(l — m)) and p = 0ifm = l — l/(3n 7 ). 
The entropy and the temperature for local equilibrium satisfy the Gibbs equation 

TdS = (pm + p M )d (J^j + jjdp M (28) 

where N is the number density, and satisfy the conservation equation 

N + 3HN = (29) 

from which we get the solution 

N(t)=N (J-^\ " (30) 

or equivalently Na (t) = const. 

Using (|28|l and l|30(l we obtain the well known evolution equation for the entropy (neglecting the heat flux and the 
shear viscosity): 

* = -^. (31) 
NT v ' 

^From eqs. ifTTjl. (THE)) . (23) , (23 and (301 we have 

(3n-l)-l/-y(l-ro) 
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For a reasonable physical behaviour we must satisfy S > 0, which implies the condition m < 1 — l/j(3n — 1). 

The total entropy S in a comoving volume is defined by S = SWa (t). Then, by considering Eqs. 1)30(1 and 11310 we 
may write the growth of total nondimensional comoving entropy over a proper time interval from tj until t / as |l lj : 

3 /**> nHa(t) , 

1 v -eft. (33) 



k B J ti T 

Now taking respectively t% and £/ as the beginning and the exit time for the inflation period of the universe, and 
from eqs. 1)16(1. 1(17(1 . 1(14(1 . 1(24(1 . ((27(1 and ((33(1 . we obtain for the increase in total nondimensional entropy in the 
comoving volume a (t) the following expression: 

3n t-t^ — r / , \ 3n— - 



The typical values for the beginning and ending times of inflation are U « 10 -35 s and tf w 10~ 32 s respectively. In 
our following numerical considerations we shall take the reference time to equal to the ending one, i.e. to = tf = 10 -32 
s. Considering that the universe at the end of the inflation period exits to the radiation era, we can constrain some 
of the constants of integration of the above formulae. Effectively, we know that the temperature of the universe at 
the beginning of the radiation era was T = 10 14 GeV — 1.16 x 10 27 K [29j. Then we have that the temperature of the 
end of inflation must be Tf = 1.16 x 10 27 K. On the other hand we know that during this radiation period p = a r T A 

7T 2 fc 4 

is valid, where a r = 5 3 ^ a = 7.56 x 10~ 15 J m~ 3 K~ 4 . From this we conclude that at the end of the inflation period 
(or at the beginning of the radiation era) we have p f as 10 93 Jm -3 . This implies that in Eq. I(14|) r = 1/4 for the 
temperature, i.e. 7 = 4/3. 

Then we can summarize these typical values for the inflation period as [Til I30I l32| : 

^wl0~ 35 s; t / «10" 32 s; Oj w ct u 
T/« 10 27 K, 7 = 4/3, p w xl0 93 J/m 3 . (35) 

The e-folding parameter Z = ln[a(t f) / a(ti)] for the power law inflation 1(16(1 takes the form 

Z = nln{^\. (36) 

It is well known that for solving the problems of the standard model in cosmology we must have Z « 60 — 70. Thus 
from (JTU), pijl. (SU and (J2H we have 



4a i e 3 ^(p{ f ) 3 / 4 
3fc B /i 



x 10 



_>(, 



(l- (10- 3 ) 3 ""^)) , (37) 



where we have considered the relation a 3 = afe 32 following from the definition of Z . We see from Eq. ((37(1 that, if the 
inequality 3^7(1 — m) >> 1 (as we have obtained above) is satisfied (in this case S/ >> Sj), we obtain the accepted 
value for the total entropy in the observable universe [23, H(j 

S« 10 88 . (38) 

Thus this model can account for the generally accepted entropy production ((38|) via causal dissipative inflation, 
without re-heating. 

Thus, we describe a dissipative accelerating period for the universe which resembles a warm inflationary scenario. 
To see this more clearly let us consider eq.(jBJ). Introducing the scalar inflaton field \ m such a way that we define 
Pm = \ jd + V(x) and pm = \x 2 ~ V(x)j where V{\) is the effective potential associated with the inflaton field, we 
obtain the following expression: 

(3H + T)x + V'(x) = 0, (39) 

where we have introduced T such that Fx = 3Htt/x- This term describes the density energy dissipated by the x ne ld 
into the thermalized bath. We have also used the slow-roll-over condition given by x <C (3 if + T)x (or equivalently, 
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p M -C +3H('jp M + 7r)/(7 — 1)). Expression (|39|l is the basic equation which describes a warm inflationary universe 
model j23|. One of the interesting characteristics of this model is does not need an intermediate reheating period as 
in usual inflationary models, since the transformation of vacuum energy into radiation energy occurs throughout the 
inflationary period and thus the inflationary epoch smoothly terminates into a radiation dominated regime 311. 

At the end of the inflationary period, when the universe enters into the radiation dominated period, we have found 
that the growth of entropy which results is £/ w 10 88 . This result gives an answer to the total entropy problem, 
whose question is why the total entropy in the observable part of the universe is so large. Certainly this problem is 
one of the various "puzzles" found in the standard Big-Bang cosmological model. 

A special case is obtained here when the parameter m = 1/2. In this case we obtain a constraint for each and all 
of the constants except p° . Effectively, substituting the expression p M = p° (io/t) 2 into Eq. 112311 we obtain 



a 



'3/ Kn(3n 7 — 2) 



9(1 - 7 /2)n 2 + 6n - 2/7 



(40) 



Now the relaxation time should satisfy the condition l|15|) . But r = ap^f and considering (|11|) we get a « y/3/n. 
This finally implies that 



n » I/7, 



(41) 



which is satisfied for a power law inflationary universe, since there n is greater than one. We could get an estimation 
of the parameter n by using the result (j38(l together with Eq. (|37Jl in which m = 1/2. We get that n > 1.6 for Z < 70. 
Case m = 1: For this value of m we have to start from Eq. (|23|l which gives the solution [2]} 



= p° e" 27t/Q 

Pm P m e 



t 



-7(3n+2) 



where n 2 = 2 7 /9 and p° is an arbitrary constant. Unfortunately this solution implies a decreasing entropy 



S(t) 



p N ( 



W 



\Vt ft 



-It I a 



const 



(42) 



(43) 



which is a non-physical behaviour. 



IV. DE SITTER FRW CAUSAL INFLATIONARY MODEL 

We also can consider the de Sitter kind of inflation models, for a closed FRW with flat dynamics. In this case 
the scale factor is given by a(t) — e Hot , where Ho = const, and the solution for the scalar field and its potential is 
obtained from Eqs.lJSJ, © and 11 21 and is given by 



4>(t) = 4>o e 



-H t 



(44) 



and 



V&) =v (f 



(45) 



respectively, where 0o = ±y/2/n / Hq and Vq = 2/n. 

Now we have H — Hq or equivalently, by Eq. (|1 1|) . p M = const. Thus from we obtain that the viscous pressure 

is given by tt — — 7 p M and then from l|31|l we see that always S > 0. Now Eq. Ijl3|l yields H — 7 p /(3a(l — 7 /2)) 
and considering Eq. I|ll|) . we have for m =/= 1/2: 



H = Hn 



3 m a(2- 7 ) 



2k 



l/(l-2m) 



(46) 



with 



(, = Co = a - 

K 



T = T Q = p[-H t 



( 7 -l)/ 7 



(47) 
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For m = 1/2 we conclude that H = H, where H is an arbitrary constant and 

27 

a = ■ . 

V3^(2 - 7) 

Here 



(48) 



(7-l)/ 7 

H = a X j- K H,T = T = n\- K HZ) (4!)) 



3 ~ „ „ /3 - 



It is easy to show that in both cases, i.e. for any value of to, we have that the relaxation time is given by 

T 3fT(2- 7 )' 

Then from Eq. (|15|) we have 



(50) 



1>\ (51) 



in accordance with the results obtained in [14( • 

Now considering that a(t) — e Ht , Z — H(tf — ti) and from l|3*3*)l we obtain 



where c has been restored. From Z ss 60 and l|35() we have that £/ >> Si and obtain that E/ w 10 88 in agreement 
with the accepted value at the end of the period of inflation. 



V. CONCLUSIONS 



In this paper we have described a closed model universe composed of a viscous matter and a dark energy components. 
The dark energy together with the curvature term "confabulate" so that they result in the dynamical Friedmann 
Equation for a flat FRW universe, as is expressed by Eq. (|llfl . On the other hand, the viscous matter component 
gives rise to an inflationary universe model at non- vanishing temperature. 

We have used the anzats (|16[) . i|14|) and l|24|l which introduce the following five parameters: n, a, to, r and /1. There 
exists a range of these parameters that gives an appropriate value for the entropy of the universe. In fact we saw that 
n > 1, r = (7 — l)/7, where 7 is the equation of state parameter, m < 1 — l/7(3n — 1) for to ^ 1 and a being an 
arbitrary constant. For a given value of a we would get a value of p° M by using Eq. I)25[l. From these ranges we have 
gotten an appropriate value for the entropy of the universe E « 10 88 . We could say the same for to = 1/2, where 
now a is given by Eq. (|4U|) . which should satisfy a « \/3/k and n >> I/7. Notice that, in this latter case, the 
matter component corresponds to quintessence, since the pressure becomes negative. Thus, the value obtained for the 
entropy in our model agrees with the expected value obtained from usual inflationary re-heating, which accounts for 
the generally accepted entropy production, via warm inflation|23|. The radiation period, after inflation, is continuously 
produced by the decay of the inflaton scalar field. In this way this field should be coupled to ordinary matter and, 
contrary to usual inflation, primordial density fluctuations are originated from thermal fluctuations rather than from 
quantum fluctuations of the inflaton [24^. 

Our results have been determined both for power-law (i.e. a ~ t", n > 1) and for de Sitter (i.e. a ~ e Hot , 
Hq =const.) inflationary universe models. 
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